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We investigate spherically symmetric spacetimes with an anisotropic fluid and discuss the existence 
and stability of a dividing shell separating expanding and collapsing regions. We resort to a 3 + 1 
splitting and obtain gauge invariant conditions relating intrinsic spacetimes quantities to properties 
of the matter source. We find that the dividing shell is defined by a generalization of the Tolman- 
Oppenheimer-Volkoff equilibrium condition. The latter establishes a balance between the pressure 
gradients, both isotropic and anisotropic, and the strength of the fields induced by the Misner- 
Sharp mass inside the separating shell and by the pressure fiuxes. This defines a local equilibrium 
condition, but conveys also a non-local character given the definition of the Misner-Sharp mass. By 
the same token, it is also a generalized thermodynamical equation of state as usually interpreted 
for the perfect fluid case, which now has the novel feature of involving both the isotropic and the 
anisotropic stress. We have cast the governing equations in terms of local, gauge invariant quantities 
which are revealing of the role played by the anisotropic pressures and inhomogeneous electric part 
of the Weyl tensor. We analyse a particular solution with dust and radiation that provides an 
illustration of our conditions. In addition, our gauge invariant formalism not only encompasses 
the cracking process from Herrera and coworkers but also reveals transparently the interplay and 
importance of the shear and of the anisotropic stresses. 
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I. INTRODUCTION 

The universe close to us is inhomogeneous exhibiting 
structures at different scales that are the result of the 
non-linear collapse of overdensities, and below certain 
scales these structures seem to be immune to the overall 
expansion of the universe. On the other hand, this pic- 
ture reveals two different dynamical behaviours that we 
wish to describe by a global general relativistic solution. 
This solution must exhibit expansion on the large scales, 
and infall at smaller scales, eventually producing bound 
structures. It is the understanding of the interplay be- 
tween collapsing and expanding regions within the theory 
of general relativity (GR) that we aim to address here. 
This issue is connected to the general problem of assess- 
ing the influence of global physics into local physics [l|, Q , 
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as well as to the approach to non-perturbative backre- 
action through model building^j, i^]. Another related 
problem is that of recollapsing 13-13 

In a previous paper fio'l, we have obtained local con- 
ditions for perfect fluid solutions to collapse within an 
otherwise cosmologically expanding background (also in 
[Til . IT^ ) . We have characterised the locally defined sep- 
arating shells between the collapsing and the expanding 
regions. 

In the present paper we wish to deepen our under- 
standing of the problem under consideration by overcom- 
ing the limits placed by the consideration of a perfect 
fluid. While such a description of the matter content is 
justified when one deals with an equilibrium configura- 
tion, the consideration of non-equilibrium states requires 
a more gen eral viewpoint where anisotropic stresses are 
present [1J|- Indeed, only for a static equilibrium does 
one expect to find identical pressures along all the spa- 
tial directions, with the exception being spatially homo- 
geneous models (in the latter case one might even ques- 
tion whether this should hold only for the isotropic case). 
So one should envisage different directional behaviours, 
which is precisely what should be expected both from 
collapsing or expanding regions within spherically sym- 
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metric models, since the radial and transverse directions 
behave differently. 

In the present work we investigate spherically sym- 
metric spacetimes with an anisotropic fluid, but no heat 
fluxes since we want to concentrate on the role of the 
stresses, shear and intrinsic curvature regarding the prob- 
lem under consideration. We leave the role of heat fluxes 
for a subsequent work. As in our previous works jiol . [T^ 
we resort to a 3 4- 1 splitting, which allows for a full met- 
ric describing both collapse and expanding regions, and 
thus avoids having to deal with the matching problem. 
We thus make use of a single coordinate patch. This non- 
perturbative approach relies on the use of the formalism 
which has been developed in a remarkable series of papers 
by Lasky and Lun using Generalised Painleve-Gullstrand 
(hereafter GPG) coordinates [l5l - [l7| . We assess the exis- 
tence and stability of a dividing shell separating expand- 
ing and collapsing regions, in a gauge invariant way. The 
local conditions that we find generalize our previous re- 
sults, and relate intrinsic spacetimes quantities to quan- 
tities characterizing the matter source. This happens 
through a generalization of the Tolman-Oppenheimer- 
Volkoff equilibrium condition, which, itself, is a gener- 
alization of the corresponding isotropic generalized TOV 
condition found in [10]. Our condition establishes a rela- 
tion between the pressure gradients, both isotropic and 
anisotropic, and the strength of the fields induced by 
the Misner- Sharp mass inside the separating shell and 
by the pressure fluxes. This defines a local equilibrium 
condition, but conveys also a non-local character given 
the definitions of the Misner-Sharp mass, and of the en- 
ergy function E (see definition in Eq. (jll.ip below). By 
the same token, it is also a generalized thermodynamical 
equation of state as usually interpreted for the perfect 
fluid case, which now has the novel feature of involving 
both the isotropic and the anisotropic stress. 

In addition, this approach has allowed us to express the 
Einstein field equations as a dynamical system involving 
scalar invariants and local quantities. This formulation 
reveals the fundamental roles of combinations of expan- 
sion with shear and two sets combining the electric Weyl 
with anisotropic stress scalars that are discussed in their 
flow evolution, relation to curvature and impact on shear 
evolution. 

To illustrate our results we analyse a particular solu- 
tion with dust and radiation. Such a solution stems from 
the work of Sussman and Pavon [IS^ where, albeit the 
generality of their initial formalism, they analysed only 
the thermodynamic aspects of the spatially fiat spherical 
solution. We find the conditions caracterising the matter 
and radiation content to fulfill the existence of a sepa- 
rating shell. In turn we also obtain the non-fiat elliptic 
solutions. 

On a different context, Herrera and co-workers (l9| 
have studied small anisotropic perturbations around 
spherically symmetric homogeneous fiuids in equilibrium. 
They concluded that this may lead to instabilities that 
result in the "cracking" of boundary surface of compact 



objects in astrophysics. We recover their results within 
our gauge invariant formalism, which not only confirms 
the important role of the shear and of the anisotropic 
stresses but also reveals transparently their interplay and 
how they trigger the cracking process. 

An outline of the paper is the following: in Section 
(|TT| the GPG formalism of Lasky and Lun and the 3 -f 1 
splitting is revised. We also define gauge invariant kine- 
matical quantities. In Section pill) we discuss the exis- 
tence of a shell separating collapse from expansion and 
give general dynamical conditions. In Section (jivp we 
present illustrations with a dust plus radiation solution 
and with the relation between the separating shell and 
cracking. Section (|V| gives a discussion of our results. 

We shall use = SttG, c = 1 and the following index 
convention: Greek indices a, 13,... = 1,2,3 while Latin 
indices a, b, ... = 0, 1, 2, 3. 



II. 3 + 1 SPLITTING AND GAUGE 
INVARIANTS KINEMATICAL QUANTITIES 

We set the basic equations in generalised Painleve- 
Gullstrand coordinates following_^the formalism devel- 
oped by Lasky and Lun (LL) [1^ [l3|, while adapting 
their derivations for our standpoint which is concerned 
with the collapse within an underlying overall expansion, 
rather than by collapse on its own. 

A. Metric and ADM splitting 

We assume that the fiow of the fiuid is character- 
ized by the timelike, normalised vector Ua — aVai — 
[— a, 0,0, 0] {nau"' = —1), defining with its lapse N — 
a and its radial shift vector N^^ = (/3,0,0), and an 
evolution of the spatially curved three-metric '^g^,^ = 

disLg (^j^ , , s'm^ 9^ . Consequently we write the 

spherically symmetric line element as 

ds^ = -a (i, Rf df + ^ ^ {(3 (t, R) dt + dRf 

+ r{t,Rfdn^, (II.l) 

which adopts the GPG coordinates of Ref. [Tt'I (dfi^ := 
d6^ -|-sin^ 9d(j)^). Notice that the areal radius r differs, in 
principle, from the R coordinate to account for additional 
degrees of freedom that are required to cope with both a 
general fluid that includes anisotropic stresses and heat 
fluxes. However, in our case we shall ignore the heat 
fluxes 



Thus, as discussed below, it becomes possible to restrict to R = r. 
However, we will restrain here from doing this identification to 
maintain the full generality in the following equations. 
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Performing a ADM 3+1 splitting [MllalSdl we use the 
projection operators along and orthogonal to the flow 



{11.2) 



where h is the 3-metric on the surface normal to 
the flow. Those projectors are also used for covariant 
derivatives: Along the flow, the proper time derivative of 
any tensor X°-^^ is 



X 



Crab 



cd 



_e \rab 
n X 



(II.3) 



and in the orthogonal 3-surface, each component is pro- 
jected with h (the overbar denotes the covariant deriva- 
tive and tensor full orthogonal projection) 



(II.4) 



Then the covariant derivative of the flow, from its pro- 
jections, is defined as 

na-b = N^Ua-c + n^.l = -ribha + ^&hab + (Tab 

+ UJab, (11.5) 

where the trace of the projection is the expansion of the 
flow, Q — — n":^^ the rate of shear aab is its sym- 
metric trace-free part and its skew-symmetric part is the 
vorticity ujab- 

On the other hand, we consider an energy-momentum 
tensor 



Tab _a b , o uo,b , -t-tab 

— pn n + F n +11 



(II.6) 



where p is the energy density, P is the pressure and n"*" 
is the anisotropic stress tensor. Ii"-^nb = and n"a = 0, 
i.e., the anisotropic stress IT^'' is orthogonal to n" and 
traceless. 

The spherical symmetry implies that all the quanti- 
ties Xap = ha°'hp°' Xab sharc the same spatial eigen- 
directions characterised by the traceless 3-tensor P"^ = 
diag [—2, 1, 1], such that 



Xap = W{t, R) Pap . 



(11.7) 



In fact, one can define in 4 dimensions the traceless pro- 
jector as Pab = hab — /tg ^^^d , admitting from the bulk 
all the properties of the 3-projector in the hypersurfaces. 
Because of spherical symmetry, we then can decompose 
any spatial 2-tensor into its trace and traceless parts: 
Xab = Kh^Xcd = ^ V + ^^ib, with its trace being 
X ~ X"^ and its traceless eigenvalue being There- 
fore, we have the following decompositions for traceless 
quantities: 



• For the anisotropic stress 

Uab^n{t,R)Pab 



(II., 



• For the shear tensor (traceless extrinsic curvature) 
aab = ait, R) Pab (11.9) 



• For the trace-free, 3-dimensional Riemann tensor, 
which measures the departures from constant spa- 
tial curvature 

^^^i?a0 - I ^'\ga/3 ''^i? = qit, R) Po.fi (11.10) 



• For the trace-free Hessian of the lapse function 
^ (d^D^ - i '^"'^g^^D^Df^ a = e(i, R) P^^ (11.11) 



• For the electric part of the Weyl tensor we have 

Eab ^ ^{t, R) Pab. (11.12) 

There is no magneticpart of the Weyl tensor due to the 
spherical symmetry [2l| , and thus the models fall into the 
class of that has been dubbed silent universes 
Another consequence of the spherical symmetry is that 
the flow is irrotational, ujab = 0. 

B. The Einstein Field Equations 

It is well known that the ADM approach separates 
the ten Einstein's Field Equations (EFE) into four con- 
straints on the hypersurfaces and six evolution equations. 
Spherical symmetry reduces them to 2+2. 

The EFE can then be written as a set of propagation 
equations: the trace and tracefre^l orthogonal contrac- 
tions of the EFE, the double orthogonal contracted and 
flow projected the Bianchi identity (once with the trace- 
free orthogonal projector read 



3_R 
2 



-2£„e + + 9a^ - -D^'D.a + 3k^P - 3A, 

a 



Cntt = — aQ + e ~ q + K H 

2 2 

Cn^ = - Y^--^ - Y (p + ^ - 2n) a 
3. + ^n) (| + « 



(11.13) 
(11.14) 



(11.15) 



They are accompanied with spacelike constraints: the 
gauge invariant radial balance, which proceeds from the 
cross projection of the EFE, and the tidal forces, obtained 
from the double orthogonal contracted, acceleration pro- 
jected Bianchi identity (again, once with the tracefree 



^ Note the sign differences in front of the Lie derivatives terms com- 
pared with [16] ; otherwise the Raychaudhuri equation restricted 
to the FLRW case does not get the usual sign for H. 

^ In other terms, contracted with h'^^n°'P%, the Bianchi identities 
yield the Weyl evolution. 
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orthogonal projector j3 yield 



and 



= —3a- 



4-77 



(p + 3n)' 



-S' - 3 S 



-n 



(11.16) 
(11.17) 



M K 



1 



Finally, the Hamiltonian constraint reads, in the presence 
of a cosmological constant, 



^^^R + -e^ - 6a^ =2k^p + 2A. 



(11.18) 



From the twice contracted Bianchi identities we also de- 
rive, along and orthogonal to the flow. 



= (i?^ + n^) (n,fc + h,kP) 
- Hi [QP + ma] . 



(11.19) 

[p-{p- 2n)] 

(11.20) 



The latter equation gives the heat fluxes evolution [T7| . 
which we set to zero here, since we are restricting our 
analysis to the case where these fluxes are absent. We 
thus have 

a' , r' 

o^-{p + p-2n) — (F-2n)+6n-. (11.21) 

a r 



The inspection of the system of equations (jll. 131111.21]) 
thus tells us that the anisotropic stress shows up in all 
but Eqs. pLT3| . (|IlT6|) . and (|TlT8| . This reveals the 
importance of the anisotropic pressures in explicitly con- 
tributing to the evolution of the shear, the electric part of 
the Weyl tensor, and of the lapse function a [19, 24-2§|. 

It is worth noticing at this point that we have included 
the cosmological constant A for the sake of completeness. 
However, none of the results that follow will depend on 
its presence. Indeed, we can without loss of generality 
make A = 0, or alternatively absorb it into p and P. 

Introducing the Misner- Sharp mass (s^l and following 



M' = ^pr'r' 



(11.22) 



it is possible to derivc[f| 
2M 



+ {\ + E) {r'y -l + -Kr^ (11.25) 



^Clr^- + -{P-2n)r^{l + E) -r' - 3 Ar . 

(11.26) 

This allows us to extend the generalization of the TOV 
function made in [lo| to the case where anisotropic 
stresses are present: 



gTOY = -Clr. 
Since, in the absence of heat fluxes we have 
1 



a 
a 



(p- 2n)' - 6n- 

r 



(p + p - 2n 

then Eqs. (|lL26l) and (llL27l) become 
M k2 

gT0Y=-Clr^- + -{P-2W)r 



(11.27) 



(11.28) 



(1 + ^) r' 
{p + P- 2U) 



[p - 2n)' - 6n- 



Ar. 

3 



(11.29) 



This tells us that, when going from the isotropic perfect 
fluid to the case of an anisotropic content in the above 
equations, we have to replace P by P — 211 and introduce 
an extra term related to anisotropic stresses. 



III. GENERAL CONDITIONS DEFINING A 
SHELL SEPARATING EXPANSION FROM 
COLLAPSE 

We now derive the generalized local conditions for the 
existence of a separating shell at r = r^,. First, we require 
a stationarity condition on the shell 

(/:„rj^ = ^ + (l + i?.) {r:f-l + ^rl=Q, [111.1] 

and second, we need an equilibrium condition to be sat- 
isfied on the shell 



-Ll = —1+ (P,_2nOr.-(l + £;.) ^ri 



a* 



A 

— Ti, 

3 



(III.2) 

Indeed, from Eq. (IIL29|) . the gTOV^ equation of 
state for the stationarity of the separating shell becomes 
now 



Or more precisely, the Bianchi identities contracted with 
/ign''PJJ, so they yield the Weyl constraint. 
^ By analogy with the perfect fluid case, it is also possible to derive 



r' Cr,E = 2{\ + E) 



(Cnr)' ■ 



C„M : 



-— r2 (P-2n) C„r. 



(11.23) 
(11.24) 



1 



(p + P- 2n) 

M 



r' 

(p - 2uy - 6n- 

r 

!^{p-2n) - iA 



<M. - 1 Ar2 

r 3 



(III.3) 



Thus the existence of a spherical shell separating an ex- 
panding outer region from an inner region collapsing to 
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the center of symmetry, depends essentially on two con- 
ditionsH The former (|IIlT|) amounts to the vanishing of 
the kinetic energy of the shell, and establishes the precise 
balance between the analogues of the total and potential 
energies at the dividing shell. The latter condition (|III.2p . 
combined with the former (|III.ip , is the generalization of 
the TOV equation for the present case, and is necessary 
for the equilibrium of the shell. There are noticeable dif- 
ferences with respect to the original problem in the form 
of the TOV equation (sil . Is^ . The isotropic pressure gra- 
dient P' is replaced by (P — 211)', the gravitational mass 
p+P \s consistently traded into (p+P— 211), and there is 
a new additional term,— 611^, involving the anisotropic 
stress n and hence reflecting its additional contribution 
to the balance of pressures and forces per unit mass. It 
is worth stressing that our result does not rely on the 
assumption of a static equilibrium of the spherical distri- 
bution of matter, and consequently does not assume that 
all the internal spherical shells are constrained to satisfy 
the TOV equation. Here the generalized TOV equation 
is just satisfied at the dividing shell. On the neighboring 
shells it won't be satisfied, and these shells will either 
be collapsing or expanding since they are not in equilib- 
riuml3 Moreover, the generalized TOV function depends 
on the spatial 3-curvature in a more general way than 
the original TOV function. 

It goes without saying that, from the conditions pil.ip 
and (IIII.2II . it is straightforward to realize that the ab- 
sence of pressure gradients between the neighboring shells 
prevents the existence of a separating shell in the spa- 
tially homogeneous FLRW models. 

Since we have 



e 



(III.4) 
(III.5) 



we see that the turning point condition piLlf) does not 
imply necessarily the vanishing of the expansion nor 
of the shear, but it rather means that these quantities 
should satisfy 9* = — Sov, at the separating shell r = 
as there 



+ 3a* =0 
=0. 



(III.6) 
(III.7) 



We emphasize again that A is written here for the sake of gener- 
aUty but is not required for the conditions to hold. 
We won't consider here the possible case where the inner shells 
move outwards and the outer shells move inwards, so that shell 
crossing results. Here we are just interested in characterizing 
the converse situation where the inner and outer shells depart. 
The occurrence of shell crossing in inhomogeneous models with 
anisotropic pressures is discussed in |33|. 



If one of O or a were to vanish at this locus we would 
then have the other quantity vanishing as well. This limit 
case corresponds to the total staticity of the separating 
shell. 



The relation in dynamics between non-local and 
local conditions 



Although the conditions pil.l[) and (jIII.3[) that char- 
acterize the separating shell hold locally, at r = r*, they 
involve non-local quantities, namely M and E. Indeed, 
from the construction of M, Eq. pi.22p . we see that the 
profile of the distribution of matter inside the separating 
shell is taken into account. 

It is however possible to find local conditions involv- 
ing local, rather than non-local quantities and this is ad- 
dressed in what follows. 

Given Eq. pil.4p it is possible to relate the condition 
pi.25p to the Hamiltonian constraint pi.lSp that gen- 
eralizes the Friedman equation. With that purpose, we 
recast the latter (also known as the Gaufi-Codazzi equa- 
tion, obtained from ^n'^inl'Gab) as 



e 



6 



A „ /e 

- + 2a — 
3 I 3 



(III.. 



so that we conclude that 

(III.9) 

In parallel, we also wish to clarify the relation between 
the gTOV function, expressed with the gauge invariant 
ofEq. (fimi) . 



gTOV 




(III.IO) 



and the "generalized" Raychaudhuri equation, obtained 
from contracting the Ricci identity with the combination 
of projectors ~\ {2h'"' + P"") ^^^ 



e 



1 

3a 



--e + --D'^Dka 



s n , 

2/3 



(P + 3P) 



(III.ll) 

It is interesting to relate '^'^^P to E from its metric ex- 
pression (jA5|) 



(1 + ^) - +1^-2 



1 ^VTTe 



VT+Er"^ 
(III.12) 



We see that the separating conditions (IIII.ll) and 
PII.3P now translate into (from Eq. IIII.8I) 



(3)P, 



(III.13) 
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and (from Eq. IIII.lip 
1 



3a 



D'^Dkalr^ ^- — {p + 3P)\ 



K 



A 

a" ' 

(III. 14) 



The former of these equations reveals that the station- 
arity condition requires ^^^R > 0, when p, A > 00 It 
no longer expUcitly involves the Misner-Sharp mass M^,, 
but just the local energy density p^, . The latter condi- 
tion emerges from the generalized Raychaudhuri equation 
pil.lip and, besides involving local quantities defined at 
r = r^f as well, it reveals that the important role of the 
pressure gradient of Eq. (IIII.3I) is now translated by the 
Hessian trace and traceless eigenvalue on the left-hand 
side of Eq. (IIILTl) . 



B. Non- locality around the shell 

It is possible to express the expansion scalar Q in terms 
of the areal radius and its Lie and radial derivatives: 



e 



and from it to derive 



iCnr)' 



(III.15) 



(III.16) 



This expression reveals that, in the inhomogeneous spher- 
ical models, the expansion scalar is not just the loga- 
rithmic derivative of the spatial volume along the timelike 
flow, unlike what happens in the spatially homogeneous 
Friedman-Lemaitre-Robertson- Walker (FLRW) models. 
Indeed, we see that it rather contains the logarithmic 
Lie derivative along the flow of the areal radius r and of 
its radial gradient r'. 

From pil.l6p . upon integration and choosing a fixed 
fiducial areal radius tq defined as tq — r(t,Ro(t)) = est, 
we obtain 



1 



6 rMr 



1 



(III.17) 



This result shows that the turning point condition at r^, 
yields 



(III.18) 



The integral on the right-hand side vanishes if the ini- 
tial parameter [r^£„r] vanishes at some interior value 
ro < r*. This requires the vanishing of the expansion 



at some intermediate value of r, tq < r < r^,. since 
it has to change signs within the interval of integration 
(we assume that no shell crossing occurs in that range). 
Differentiating equation (|III.17p with respect to the fiow, 
we obtain 



2Cnr 



J ra / 



[Cnr] e - - [Cnr] 



gTOV, 



(III.19) 



where r denotes proper time. This is the equation that 
generalizes the Eq. (3.27) of llOll and that corresponds 
to Eq. (21) of di Prisco et al. [3J|. It corroborates once 
again the claim of a non-locality of the radial accelera- 
tion. From Eq. (jlll.l7P we realise that this non-locality 
is inherent in the radial expansion, and is already present 
in the energy condition defining Eqs. pil.llIII.6|) and 
in our gTOV condition Eqs. pi.27lll.26IIII.2|) . since both 
implicate M which is an integral between and r-^. 

From the previous equations piI.17|) and piI.lQp we 
see that at the separating shell we have 



56 2 , 

— — r dr 

OT 



(III.20) 



Strictly speaking, when k p* -I- A > 0. 



which means that the integral on the right-hand side van- 
ishes if the term — [r^£„r]^ vanishes at an interior 
value To < r*. This shows that the vanishing of the 
proper time derivative of the expansion >C„0 occurs then 
at some intermediate value between ro and r*. In the 
case when £„ [r^>C„r] ^.^ = at the center, we recover 
the result of Di Prisco et al. [34], establishing the van- 
ishing of the radial acceleration, i.e. CnQ ~ 0, at some 
< r < Tvf. However this results is derived here in a 
non perturbative way, and in a more general way than in 
Ref. Isl. 



C. Dynamics around the shell 

In this section, we will address the dynamics of the 
system under consideration, adding various restrictions 
of interest for the rest of the paper and, in each case, 
examining the dynamics of the matter-trapped shell. 



1, Dynamical system of the imperfect fluid 

a. Governing equations The dynamical system of 
partial differential equations (PDE's) that results from 
the 3+1 splitting and the use of the local kinematical and 
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geometric quantities is given by Eqs. pil.llllll.141lll.l31 
UlLSl IIL2T1 UlTel [TlT7| and a constraint (Eq. IIII.24P on 
the Weyl tensor induced by the differences in the shear 
equation obtained by projections both from the Einstein 
field equations (Eq. III. 141 from ^^Gcd) and from the 
Ricci identities f Eg. 1111.2^ from the projection— |P°^n''). 
We restate the whole system as 



1 

3a 



(p + 3(p-2n)) 



e 



£„a = — a8 + a ( "^ + '^ 



A 



(III.21) 



(III.22) 



what follows. E + -^n acts as a source for density inho- 
mogeneities as seen in Eq. (jIII.28|) . From Eq. pil.24p we 
see that this is related to the 3-curvature distortion of the 
hypersurfaces as well as to the distortion of the extrinsic 
curvature, as expected. The role of the other combina- 
tion is clearly revealed in the shearfree subsection IIII C 21 
that follows. 

Alternatively, one can present Eq. pi.l6p in a form 
parallel to that of Eq. pil.28p 



3 



=0. 



(III.30) 



h. Dynamics of the shell On the separating shell, 
the dynamics can be expressed from the EFE and Bianchi 
identities. It takes the form of the residual constraint 
from the Raychaudhuri equation (Eq. IILT31 + IILT8| /21 /6 
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3PJ- 
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III.31) 
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(III.23) 
(III.24) 



- + 2a - 



r' a' 
{P - 2X1)' =6n (p + P - 2n) — , 



3a- 



((' 



-n 



(III.25) 
(III.26) 

(III.27) 
(III.28) 



In this formulation the equations revea0 the fundamen- 
tal role played by some combinations of gauge invariant 
quantities like expansion and shear, electric Weyl and 
anisotropic stress. In the latter case, they emerge in two 
different combinations that play different and important 
roles in the governing equations, as we will discuss in 



and the "generalized" Raychaudhuri Eq. pil.lip 



3ai 6 



3PJ 



(III.32) 



The Hamiltonian constraint yield the local curvature of 
the shell, ^R^, — 2K^pi, + 2A, the momentum constraint 
governs the expansion and shear transfer across the shell, 

+ a)^ = — 3a*^, the Weyl constraint from the shear 
equations links it directly to the 3-curvature residual 



(III.33) 



the density remains conserved by Eq. (111.19^ which, with 
Eq. plLell . reads now 



e.(p + p-2n) 



(III.34) 



The Eq. pi.2ip gives a part of the gTOV staticity condi- 
tion. The Weyl constraint Eq. pi.lTp governs the balance 
of anisotropic stress and energy density across the shell. 
But most interestingly, the evolution of the electric part 
of the Weyl tensor is bound to that of the anisotropic 
stress by Ea. pil.23l) which reduces here to 



-n 



6 



(p + p-2n)^e^. (III.35) 



This is to be related with the studies on cracking by 
Herrera et al. [1£, ^] . We now restrict to shear 

free flows. 



9 Notice that Eq. IIIII.21I I can also be noted 



fcO -I- e - ( ^ + a 



^2 , , A 

- (p + 3P) + - 



(III.29) 



2. Dynamical system restricted to shear free flows 

We set out to restrict to shear free fiows as they con- 
stitute an important subcase in many studies, i.e. as in 
[2^ . [35l - [37t or even in cosmological FLRW models. 
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a. Governing equations The Eq. (III.14P reveals the 
necessary and sufficient condition for a shear-free flow to 
be 



(III.36) 



Here the combination of electric Weyl and anisotropic 
stress that govern the shear evolution appears clearly in 
its role. This generalizes the result of [23| to the case of 
non- vanishing acceleration. 

The remaining equations of the system (assuming 
IIIL36P that differ from the general case reduce to 



£„e ^-D'^D^a - |- - ^ (p + 3P) + A 



a 



£„ ( E + -n ) = - 



2 I S+yH 



2 



(III.37) 

e 

3"' 
(III.38) 
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3 

6' =0 



K p 



A, 



(III.39) 

(III.40) 
(III.41) 



Notice that Eq. (|III.39p shows that in the shearfree 

2 

case, the E + ^11 combination of electric Weyl and 
anisotropic stress relates only to the 3-curvature distor- 
tion of the hypersurfaces. It also is a generalization of 
the constraint k^II = g = 21] found in (H-fl^. More- 
over, Eq. pil.38p that governs its evolution can be re- 
expressed, using Eqs. piL39p and (|III.36I) . as 



[2g + e]|, 



so q is damped by ^ . Therefore the sign of 8 determines 
the increase or decrease of the 3-curvature distortion. Ex- 
pansion dampens the distortion while collapse enhances 
it. More importantly Eq. (jIII.4ip implies that Q does 
not depend on R, and therefore 



(9B 
at 



92 , 

1 ip 

3 2^^ 



3P)- A 



of the system reduce at that locus to 
99 

a^CnQ* - =D^Dka, - a,— {p + 3P), + a,A 



(III.42) 

where x is a function of just the time coordinate t. The 
Hessian trace is thus determined by the Friedman accel- 
eration sources and a time dependent term: ^D'^D^a = 

^ ip + 3P)-A + q^ + m, 

b. Dynamics of shear free limiting shell Since at 
we further have 9* = 0, the remaining changed equations 



= -3a^£„a^ = 



=0 

e' =0 



(HI.43) 
(HI.44) 
(HI.45) 



From cq. (IHI.44P we reahze then that ( E + -^^H 



is a constant of the motion along the flow n"' (timelike 
vector fields). 

In the shear-free case the expansion scalar throughout 
is only a function of time, and the relation between the 
local values of the electric part of the Weyl tensor and of 
the anisotropic stress does not change along the orbits of 
the shells. It is also worth noticing that the 3-curvature 
of the dividing shell is completely determined by the local 
energy density. 

A "limit" case is the case of a static initial configura- 
tion 9 = in addition to the vanishing of the shear. We 
will consider this case in the subsection on the cracking 
phenomena. 



3. Dynamical system restricted to geodesic flow 

The following case of geodesic flow is defined by no 
acceleration. The perfect fluid solutions are dealt with 
in [38] 13 This implies a' — 0, and therefore e = 0, 
D^Dka — 0, and it is advisable to set a{t) — 1, and 
use gLTB coordinates (generalisation from the Lemaitre- 
Tolman-Bondi, hereafter LTB, coordinates). As a re- 
mark, the more restrictive geodesic, shear-free flows are 
subject to 



-H 



(III.46) 



recovering the Mimoso and Crawford result [25] , and the 
subsequent discussion of Coley and McManus [2j3, . 

a. Governing equations The equations for geodesic 
flows that differ from the general case now reduce to 



r 



9 



a9 



(P + 3F) + 



A 

+ 3' 
III.47) 



-H 



(III.48) 



There, spherically symmetric, inhomogeneous models are stud- 
ied, and it is shown that there are no exact, geodesic, non-static 
perfect fluid solutions with nonzero shear. 
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(P - 2n)' =6n 



(III.49) 



This case is interesting as it corresponds to the generali- 
sation of the classic LTB model [31, 39] as well as to the 
Sussman and Pavon [I^ |40| example we will use later. 
The major difFerence from the general case is the absence 
of Hessian trace, —D°-Daa, and traceless Hessian, e, in 
the equations pH?47l) . pil.48p and (|HI.49|) . In partic- 
ular Eq. pil.49|) displays a completely different radial 
constraint: not only the inertial mass is no longer in- 
volved, as the acceleration vanishes, but also in this way 
the anisotropic stress is the only source for the inhomo- 
geneity of the pressure. For a perfect fluid the pressure 
should be spatially homogeneous, as found in [38]. 

b. Dynamics of geodesic limiting shells Further re- 
stricting to the shell r^, we have the remaining changed 
equations 

Cni^+a] =0 = -|^(p + 3(P-2n)) + (S + yH) 



A 

+ 3' 

£„a* =— S — — H 

3 V 2 

(P-2n)l =6H,^, 



(HI.50) 
(HI.51) 
(HI.52) 
(ni.53) 



The definition of the matter-trapped shell then implies 
Eq. (|HI.50P which is the local version of the gTOV, i.e., 
the local gRAY = equation. If in addition we have 
the shear free condition (IIII.46|) . we see that the matter- 
trapped shell imposes that (p + 3P) be locally constant 
or vanishing (if A = 0). On the other hand Eq. pil.53|) 
shows that the value of -I- a) in the neighborhood of 



> 



the separating shell is non vanishing and that (y -I- a 
provided 0*r^ > 0. 

c. Geodesic Misner-Sharp mass and electric Weyl 
For the geodesic fiow, we have, from Eqs. pil.2[) and 
(IIII.lip . the latter in the form of Eq. pil.47|) . a relation 
between the Misner-Sharp mass and the electric Weyl: 



M ( n\ 

y[p + 3n]-E 



(III.54) 



because of its definition (IIII.6I) , the expansion of the out- 
side region does not precisely cover the usual expansion 
region: the separating shell itself can have non-zero ex- 
pansion and thus one of the said collapsing or expanding 
region may contain the = shell. However, the choice 
was not laid on such shell because the present definition 
yields the staticity condition on that surface (Eq. IIII.3P 
which is not, in general the case for turnaround shells 
(0 = 0). 

The meaning of expansion in the terms of Sec. IIIII 
is linked with the areal radius: the luminosity distance 
of a shell to the centre. Thus the static shell keeps its 
luminosity distance to the centre while expanding regions 
appear so in the luminosity distance space. 

Isolating the Ricci curvature of spatial hypersurfaces 
in Eqs. pi.lSP and pil.8|) . however, reveals that 
both = and Eq. pil.6p require '^^i? > 0, placing the 
respective surfaces both in the positively curved region of 
spacetime, where the region between the two must lie. In 
models with negatively curved regions, those expanding 
shells will therefore be contained in the expansion regions 
defined for both expansion scalar and areal radius. The 
flat and closed background can still present expansion 
infinities in both senses, as seen in [ij], although the 
general treatment can be more complex. 



IV. APPLICATIONS OF OUR RESULTS 

A. Sussman-Pavon exact solution: radiation and 
matter 

To illustrate our results we turn our attention to an ex- 
act solution derived by Sussman and Pavon for a spher- 
ically symmetric model with a matter content consist- 
ing of a combination of dust and radiation that exhibits 
anisotropic stress, but no heat fluxes and A = [18i]. 

In order to do that we need to translate the metric 



(III.ip in the LTB form used in [18|, following their as- 
sumption of comoving, i.e. geodesic, flow (here the ve- 
locity of light c is reintroduced) 



That judiciously imposed condition translates into a flow 
without acceleration, which leads in terms of metric com- 
ponents to a = a{t). Thus, the time function can always 
be rescaled to absorb the lapse at — c. 



1. Coordinate transform from GPG to LTB 



D. Separation and expansion 

In cosmology, the expansion of the background uni- 
verse is understood as the condition on the universal fluid 
flow of O > 0. In Sec. Illli we have extended the deflni- 
tion of [l^ for matter-trapped surfaces separating expan- 
sion from collapse, however it should be explicited that, 



Canceling the metric crossed term gives a similar rela- 
tion than the perfect fluid generalised LTB formulation 
of [16], while the radial term imposes 

{pi + R^R' = Q, (IV. 1) 

i?'2=r'2, (IV. 2) 
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(with R = RcpG, the ' and ' denoting derivatives in the 
gLTB frame). Compared with the absence of heat 
flux suggests that the extra degree of freedom provided by 
a spacetime dependent lapse in the GPG frame becomes 
superfluous. The new areal radius, from Eq. pV.2p is 
reset to the GPG radial coordinate r = Rq pa ~ TZ (here 
for convenience we change notation for the GPG radial 
coordinate), yet is still spacetime dependent. Then, as 
in the perfect fluid case in [l^ , the coordinate transform 
(IrTTI) is such that /Sdt + dTZ cx dR. Taking t{T) = cj ^ 
and r(T, i?), we have then the condition 



(3dTt + OtTZ = 0, 
which becomes (in GPG coordinates) 

^ 7^ 



(IV.3) 



(IV.4) 



Moreover, Eq. pi.23p implies that, in the new gLTB 
coordinates, E = E{R). Consequently, the line element 
(III. Ill can be rewritten as 



ds'^-C^dT^ ' ^^""^^ -^^2 , .2^o2 



1 + E{R) 



dR^+r^dVL\ (IV.5) 



as in hj. 



2. Restricted dynamical equations 

The crucial anstaz adopted by Sussman and Pavon was 
the assumption that the flow is geodesic, keeping as close 
as possible to the case where dust is the only component 
present, i.e., as in the original LTB case. The Bianchi 
contracted identity pi.20p . together with the geodesic 
condition a' = <^ rt" = 0, imply that 

D'' {h,bP + n,b) - [OP + 6Ua] = 0, (IV.6) 

that is, in gLTB coordinates. 



(P-2n)'-6n- =0, 

r 



(IV.7) 



where the prime stands for differentiation with respect 
to the geodesic R. For practical purposes this amount to 
have 



and 



Mdust = M{R) , 
W{R)r,iR) 



Mrad = 



2r(T, R) ' 
so that Ea. (|n.25p . including c, now reads 



2— + + E 



(IV.: 



(IV.9) 



(IV. 10) 



From Eq. (jlV.lOp . one is led to the following so- 
lutiorP^. generalised from [l^ to encompass the cases 
where ^ 0, 



± c / At = 



M , ( E^I'^^jEr'^ + 2Mr + Wr, + Er + M 
In 



£^3/2 



E 



(IV. 11) 



Existence of a separating shell in the generalised 
Sussman- Pavon solutions 



We see that the vanishing of the right-hand side of Eq. 
pV.lOP provides one of the conditions for the dividing 
shell, while the other condition that corresponds to the 
gTOV equation will be derived from the radial accelera- 
tion 



r 

72 



M 



Wn 



(IV.12) 



or directly from combining Eqs. pi.l3[ /64 lII.14p in the 
gLTB frame with Eq. pV.lOp . So we find that 



2r3 



-iP-2U)r. 



(IV.13) 



The form of the gTOV condition in the gLTB frame (Eqs. 
III. 271 III. 261 with conditions a' — 0) can then be recog- 
nisedin Eq. pV.12p using the Raychaudhuri constraint 
(llVl3l) . 

The existence of a matter-trapped shell in the solution 
inspired by |18] requires both Eqs. pV.lOP and pV.12p 
to be zero on some r ~ r^,. This implies the existence of 

and, from Eq. pV.12p . and then Eq. pV.13p . 

W, = - Nh— ^ WiR,) < ^P,< 2U,. 

(IV. 14) 

This latter condition shows that, in order to allow the 
separating shell to exist locally, the model of Sussman 
and Pavon [18] must contain regions where the transverse 
pressures balance the radial pressure. This can be under- 
stood with Eq. pV.13p xr^: the radiation Misner-Sharp 
mass corresponds to the flux of the pressures across the 
shell. Setting Eq. pV.lOp to zero implies, using again 
Eq. (iRTTa . 



2 =-2^-K^P.-2U^)rl 



(IV. 15) 



1^ The elliptic integral, from Eq. HIV. 101 . reads ±c J dt 
r rdr _ r dY M r dX 



y/Er^+2Mr + Wri 



2^y/Y+^ £3/2 J y/x^+cst 



I- 
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2M Wr, 



Mr- 



E 



E, 



Hm Wr; 



E< 



Figure 1. Dynamical analysis of a local W < shell. The 
dynamic for a given shell (fixed M, W and E without shell 
crossing) obeys Eq. (|IV.10|I . It then behaves as a one dimen- 
sional particle in an effective potential, following We 
draw a qualitative energy diagram to illustrate the definition 
of the critical curvature/energy Eum, when it lies in a region 
of W < 0. The various cases of i?> > Eum, i?< < Eiim&nA 
E — Eiirn yield unbound, bound and marginally bound be- 
haviours. 



while Eq. (IIV.12P gives 



(IV. 16) 



so with pV.13p . the energy/curvature parameter E reads 



(IV. 17) 



Again, as in the perfect fluid case [1Q|, the separating 
shell only exists in elliptic regions {E < 0). Finally with 
Eq. PV.14|1 we have 



__W 



and thus 



E, = 



(IV.18) 



(IV. 19) 



For outward initial flows, this requires Vi < r^, thus the 
additional condition W < —M < 0. 




Figure 2. Global, qualitative, analysis yielding the separating 
shell at the intersection of E with Enm- Following we 
construct initial cosmological conditions with a.W <Q region 
such that Eqs. (|IV.10[I and (|IV.12[I are simultaneously, that 
is _E = Eiim and W — —M, so the intersection of the curves 
gives a global dynamical separation. 



velocity flow, FLRW outer behaviour (M, W ^ r^, 

r— f oo 

E ^ r^) and an intermediate W < —M < {) re- 

r— )-oo 

gion can be qualitatively obtained (Fig. (2). Then using 
Eiim — ■^Tj and choosing initial velocities such as E 
crosses Eum in the < region, one gets the global 
separation (see Fig. Although we allow for a region 
where the Misner-Sharp mass of the radiation fluid is 
negative, we remind that only the total density of the 
fluid is actually meaningful and point out that we should 
keep 



M' ■ 



+ Pr), > 0, BO M + W >0. 



This implies that only a static global separation can fulfill 
both W < -M and M + W >0: W = -M, obtained by 
crossing E with Eum at the radius where M = —W, as 
shown on Fig. ^ Note that the energy conditions does 
prevent initial conditions with inward going initial flow 
in the neighbourhood of the global separation if the no 
shell crossing condition is to be maintained. In that case 
it is allowed to have initial radius outside < but 
then the shells just outside the separating one should be 
ingoing and unbound. This would result in shell crossing 
after some time in a symmetric way as found in (l4| in 
the case of the analysis of a ALTB model. 



4- Dynamical analysis and global shell 

As for the examples of Ref . [l3| , a dynamical analysis 
of Eq. (jIV.lOp can be performed in the regions where 
W <0 (see Fig. [1]), required by Eq. (IIV.14p . Initial con- 
ditions with cosmological outwards initial areal radius 



B. Cracking phenomenon of Herrera and coworkers 

We find a second illustrative example of our results in 
the concept of cracking put forward by Herrera and col- 
laborators [41j:] whereby a static spherical configuration 
is unstable to anisotropic perturbations and "cracks". 
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In order to discuss this concept within our framework 
we have to consider the set of EFEs as in Sec. IIII C li s 
Eqs. (|III.2ip to (|III.28|) . 

From these equations we see that the shear plays a cen- 
tral role. From the shear propagation equation (|III.22p 
we realize that if the shear were to vanish initially, any 
deviations from constant curvature given by the term 

2 

£—(!] — ^n) would indeed make the shear become non- 
vanishing at any later instant. 

To recover the cracking phenomena envisaged by 
Herrera and collaborators we start assuming a static, 
isotropic, shear-free initial configuration. Thus we put 
= 0, n = and a = in some region of the initial hy- 
persurface where we assume also that p and P only vary 
slowly. Then we can assess any future deviation from 
that configuration using the restriction to this hypersur- 
face of Eqs. (|III.2ip to (|III.28p . The governing initial 
equations, for A = 0, 9 = 0, 11 = and a = 0, take the 
form 

CnQ=-D^Dka-'^{p + W), (IV.20) 
a 2 

£„a=e-E, (IV.21) 

(IV.22) 
(IV.23) 
(IV.24) 

(IV.25) 

(IV.26) 

. (IV.27) 



£„ ( S + yH ) =0, 



(p-2n)'^-(p + p)- 

a 



| + =0, 



The form of the Eq. pV.27|) uses the constraint pV.231 
actually 1111.24^ while the Raychaudhuri Eq. pV.20|) 
comes from the restriction of Eas. (IIII.2HIIII.22p to the 
initial configuration. Further using Ea. (|IV.22p with the 
constraint pil.24p in the derivative, one can deduce the 
relations between the values and proper time evolutions 
of the electric Weyl scalar, the anisotropic stress and 
traceless hypersurface curvature, as well as traceless Hes- 
sian scalar, shear and expansion on the initial hypersur- 
face 



a =e = n = 0, 



(IV.28) 
(IV.29) 
(IV.30) 
(IV.31) 



shear Eq. pil.22p on the hypersurface 

Cla =Cn (e + K^n) = Cn (e - 2S) . (IV.32) 

From all this, the following can be deduced: (i) the per- 
fect fluid source and Hessian combination in Ea. pV.20p 
drives, in general, the expansion away from 0; (ii) com- 
bining Eqs. PV.21|) with pV.23|) . the shear, in general, 
is also driven away from by the difference of traceless 
Hessian and curvature of the hypersurface, q, the latter 
mirroring the Weyl curvature, unless they are set equal, 
thus implying an imposed shear-free flow; (iii) anisotropy 
as well is driven away from 0, in parallel with the evo- 
lution both of the electric part of the Weyl, S, and of 
the difference of traceless Hessian and hypersurface cur- 
vature, q, as we have 



(IV.33) 



except if the flow is restricted to being shear-free and 
geodesic; (iv) the separation scalar defined in Eq. pH.6[) 
is on the initial hypersurface but will be driven away 
by 



r 



e 



(P + 3P) 



Therefore, if a shell where g = — ^ (p + 3P) exists in the 
considered region, it satisfies locally the TOV equilibrium 
condition, where forces balance, and is thus surrounded 
by shells experiencing nonzero forces. That shell satis- 
fies the conditions (IHI.6I) and pH.7[) of a dividing shell. 
As a consequence, we realize that gTOV becomes non 
vanishing in the neighborhood of the dividing shell, in- 
ducing the appearance of the radial force responsible for 
the cracking phenomena under the following conditions: 
for G*r^ > 0, the radial balance of the separation scalar 
will drive neighbouring shells to the cracking condition of 
outer shells and inner shells experiencing positive, resp. 
negative areal expansion, as seen in Eq. pi.l6[) . at some 
later time from those initial conditions. In conclusion, 
the cracking shell is a kind of separating shell. 

We can further extend our interpretation of cracking 
in this framework by considering shear-free flows. Then, 
(i) expansion is still driven away from by ■^D'^DkCt — 

^{p + 3P); (ii) anisotropy is still driven away from 
by the traceless Hessian; (iii) as seen in Sec. IIII C 21 
the shearfree condition entails from Eq. pil.4ip that 
the expansion should be of uniform sign in all spacetime. 
The departure from initial vanishing expansion by Eq. 
(IIV.20[) will give its definite sign and thus the expansion 
is always either all collapsing or all expanding, as found 



where, in addition, the traceless Hessian scalar proper 
time evolution can be obtained with the derivative of the 



An important point to be emphasized at this stage 
is that our analysis draws on the full set on non-linear 
equations and is therefore more general than that of the 
original works of Herrera an collaborators. Alternatively, 
a perturbative gauge invariant treatment of this issue can 
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be done using the formalism develo ped by [Z^l and sub- 
sequently explored by others (iol . I43l - l46l | . We will address 
it elsewhere. 

We conclude this section emphasising the importance 
of both the shear and the anisotropic stresses not only for 
the existence of a dividing shell, but also for the cracking 
phenomena of Herrera and collaborators. In the latter 
case, this confirms their claims in an alternative way. 



V. SUMMARY AND DISCUSSION 

In the present work we have considered spherically 
symmetric, inhomogeneous universes with anisotropic 
stresses in order to investigate the existence and stabil- 
ity of a dividing shell separating expanding and collaps- 
ing regions. With this endeavour we have gone one step 
further than in a previous work by considering a more 
realistic scenario where the matter is no longer a perfect 
fluid. 

This shows up to be quite important in characterizing 
the contrasting dynamical behaviours of separate regions. 
This is relevant in relation with the present understand- 
ing of structure formation as the outcome of gravitational 
collapse of overdense patches within an overall expand- 
ing universe, since there is an underlying expectation that 
the two disparate behaviours decouple. This issue, is also 
related to the assessment of the influence of global physics 
on local physics. 

In the present work we have addressed this issue 
by resorting to an ADM 3+1 splitting, utilising the 
so-called Generalized-Painleve-Gullstrand coordinates as 
developed in Refs. [l^ 13 • This enables us to follow a 
non-perturbative approach and to avoid having to con- 
sider the matching of the two regions with the contrasting 
behaviours. We have found local conditions characteris- 
ing the existence of a dividing shell which generalises our 
previous conditions for perfect fluids fl^ . One is a condi- 
tion establishing the precise balance between two energy 
quantities that are the analogues of the total and poten- 
tial energies at the dividing shell. (This amounts to the 
vanishing of the kinetic energy of the shell.) The second 
condition establishes that a generalized TOV equation 
is satisfied on that shell, and hence that this shell is in 
equilibrium, but one which now involves explicitly the 
anisotropic stress. Moreover, the former condition also 
implies that there is no matter transfer across the sepa- 
rating shell, and hence we may call the region enclosed by 
the latter a trapped matter region. The trapped matter 
is not in static equilibrium in contrast to the situations 
where the TOV equation is satisfied in the whole of the 
trapped region, and which are meant to describe stars. 

We have also related these conditions to a gauge in- 
variant definition of the properties of the dividing shell. 
These require the vanishing of a combination of the ex- 
pansion scalar and of the shear, on the shell. Finally, if 
we demand that the dividing shell is static, in order to 
define an extreme case of reference, we obtain an addi- 



tional equation of state that relates the gauge-invariant 
Hessian trace of the model with the quantity p + 3P in- 
volved in the strong energy condition, on the dividing 
shell. Naturally in this limit case, the expansion and 
shear will vanish on the shell. 

The approach followed in this paper has allowed us to 
translate the Einstein field equations in terms of non- 
local quantities, such as the Misner-Sharp mass M and 
the energy /curvature parameter as well as into equa- 
tions involving local quantities. The latter are conve- 
nient to describe the evolving behaviours separated by 
the dividing shell. This procedure led us to relate the 
energy equation to the generalized Friedmann equation, 
and likewise we relate the generalized TOV equation with 
the generalized Raychaudhuri equation. Moreover, we 
present both the equations governing the flow behaviour 
of the remaining quantities such as the shear, the electric 
part of the Weyl tensor, as well as the constraint equa- 
tions that hold for them and for their radial gradients. 
We also give the constraints and evolution applied on 
the expansion and shear combination. This allowed us 
to discuss the dynamical behaviour in the neighborhood 
of the separating shell. In particular we have obtained 
the condition for a shear-free flow that generalizes previ- 
ous results [Ullllll. 

We have considered two illustrations of our results, 
namely we have analysed the existence of a separating 
shell in the class of matter and radiation solutions put 
forward by Sussman and Pavon. We showed that, in this 
case, the existence of a dividing shell requires that the ra- 
diation exerts a repulsive role. And we have shown that 
our results allow a discussion of the emergence of the 
cracking phenomena put forward by Herrera and collab- 
orators. We have described cracking initial conditions, 
their dynamics, and showed, within our gauge invariant 
formalism, how shear and anisotropic stress trigger the 
phenomenon of cracking. Our approach also opens win- 
dows on the behaviours of the electric part of the Weyl, 
the quantities characterising the 3-curvature. We also 
recover the properties discussed in [SS^ in shearfree flows. 

In this paper we didn't do a thorough discussion of all 
dynamical possibilities offered by the system we discov- 
ered. This opens many future work. 
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Appendix A: Metric ADM scalar functions 
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For clarity, we present the scalar gauge invariants in- 
volved in the ADM formulation in terms of the GPG 
metric functions, starting with the perfect fluid terms 
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which, combined with Eq. pi.23[) . yield 



(Al) 



(A2) 



(A3) 



and from the shear evolution comparing that from the 
EFE, Ea. (jII.14p . and that obtained from the Ricci iden- 
tities, we get 
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